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Abstract. We present a discrete model of resonant scattering of waves by an open pe- 
riodic waveguide. The model elucidates a phenomenon common in electromagnetics, in 
which the interaction of plane waves with embedded guided modes of the waveguide causes 
sharp transmission anomalies and field amplification. The ambient space is modeled by 
a planar lattice and the waveguide by a linear periodic lattice coupled to the planar one 
along a line. We show the existence of standing and traveling guided modes and analyze 
a tangent bifurcation, in which resonance is initiated at a critical coupling strength where 
a guided mode appears, beginning with a single standing wave and splitting into a pair of 
waves traveling in opposing directions. Complex perturbation analysis of the scattering 
problem in the complex frequency and wavenumber domain reveals the complex structure 
of the transmission coefficient at resonance. 
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1 Motivation 

When a periodic waveguide is in contact with an ambient space, the interaction between modes of 
the guide and radiation originating from the ambient space outside the guide results in interesting 
resonant behavior. The resonance is manifest by pronounced amplitude enhancement of fields in 
the waveguide and sharp anomalies in the graph of transmitted energy versus frequency near the 
frequency of the guided mode [31 El El e.g.]. Examples abound in the physics and engineering 
literature because of the importance of these anomalies in applications to photo-electronic devices. 
Because of the exchange of energy between the waveguide and the surrounding space, we call the 
waveguide open. 

In this work, we analyze a discrete model of resonance in open periodic waveguides. The 
ambient space is modeled by a uniform two-dimensional lattice, and the waveguide is modeled by a 
periodic one-dimensional lattice, coupled to the two-dimensional one along a line. With period two, 
we consider this to be the simplest model that exhibits the essential features of the open lossless 
waveguide in air. It is not unlike the Anderson model, in which a single chain of beads interacts 
with a resonator attached to one of the beads |XQ^ [Tlj I12j . But, unlike the Anderson model, a 
periodic model, with period at least 2, admits both propagating and evanescent Fourier harmonics 
simultaneously, and this is precisely the feature of open periodic waveguides that allows embedded 
guided modes and their resonant interaction with incident radiation. 

The discrete model is useful in that it exhibits important resonant phenomena of continuous 
open waveguides while its simplicity permits explicit calculations and proofs. In particular, one 
can prove that the resonant peaks and dips in the transmitted energy reach exactly 100% and 0%, 
a phenomenon that is often observed in scattering of electromagnetic waves by open waveguides. 
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Moreover, explicit formulas illuminate the connection between structural parameters of the waveg- 
uide and the properties of the anomalies, such as central frequency and width, both of which are 
important in applications of lasers and LEDs [2] . 

In addition, we analyze a tangent bifurcation of resonances, in which resonance is initiated at a 
critical coupling strength, beginning with a single standing wave and splitting into a pair of waves 
traveling in opposing directions. 

The kind of resonance we are describing here is akin to those that go by the names of Feshbach 
resonance, Breit-Wigner resonance, or Fano resonance in quantum mechanics. The unifying idea is 
that, when one perturbs a system that admits a bound state whose frequency is embedded in the 
continuous spectrum, the eigenvalue dissolves as a result of the coupling of the bound state to the 
extended states corresponding to the frequencies of the continuum. This coupling is the cause of 
sharp features in observed scattering data near the bound state frequency [151 §XXIL6], [H [22] . 

A similar type of resonance of classical fields (as electromagnetic or acoustic) results from the 
interaction of guided modes of a periodic waveguide and plane waves originating from outside the 
guide. For this interaction to take place, the waveguide must be open, that is, it must be in 
contact with the ambient space, such as in the case of a photonic crystal slab in air. Because of the 
periodicity of the waveguide, guided modes may couple to the Rayleigh-Bloch diffracted waves and 
become "leaky", or "quasiguided" modes, or "guided resonances" [3l[Tll[2T]. These leaky modes 
are associated with the sharp anomalies in the graph of the transmitted energy across the slab as 
a function of frequency, as we have mentioned. 

Under certain conditions, a lossless open periodic waveguide can actually support a true guided 
mode — one that is exponentially confined to the guide. This can occur in one of two ways: (1) 
if, for a given frequency and Bloch wavevector parallel to the guide, the expansion of the fields 
in spatial Fourier harmonics parallel to the guide admits no harmonics that propagate away from 
the guide (this is the region below the light cone in the first Brillouin zone); or (2) there are 
Fourier harmonics that propagate away from the waveguide (the Rayleigh-Bloch diffracted waves, 
or propagating diffractive orders) but the structure admits a field for which the coefficients of these 
harmonics happen to vanish. The latter occurs, for example, in a waveguide that is symmetric 
about a plane perpendicular to it at a wavevector parallel to the plane of symmetry [IJ [20l [21] . 

We will be concerned with the latter type of guided mode, for only these can interact with 
incident radiation. Such modes are typically nonrobust and are therefore excited by small deviations 
of the angle of incidence (associated with the Bloch wavevector) or perturbations of the structure. 
Their frequencies can be viewed as embedded eigenvalues in a pseudoperiodic scattering problem 
for a fixed Bloch wavevector. Perturbation of the wavevector or the structure itself destroys the 
guided mode and the associated eigenvalue. In [19], Shipman and Venakides derived, for the 
two-dimensional case, an asymptotic formula for the transmitted field as a function of wavevector 
and frequency based on complex perturbation analysis of the scattering problem about the guided 
mode parameters. The formula is rigorous and subsumes that derived by Fano in the context of 
quantum mechanics [3]; it plays a major role in the analysis of resonance in this paper. An in-depth 
discussion of resonance near nonrobust guided modes can be found in |16] , and the role of structural 
asymmetry on the detuning of resonance is analyzed in a discrete model in |17j . 

The exposition proceeds as follows. 

§2. The discrete model. First, we describe the ambient and waveguide components of the 
discrete model and how they are coupled. We then identify the minimal space of motions in the 
ambient 2D lattice (the "reconstructible" part), which, together with the ID waveguide, form a 
closed system, decoupled from a complementary "frozen" space of motions of the 2D lattice. 
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§3. The scattering problem. The problem of scattering of plane waves in the 2D lattice 
by the ID waveguide is posed. We describe the Fourier decomposition of fields and resolution 
of scattered fields into their diffractive orders and prove existence of solutions, including in the 
presence of a guided mode. 

§4. Guided modes. We prove the existence or nonexistence of true guided modes under certain 
conditions and describe the dispersion relation for generalized (leaky) modes relating complex 
frequency to complex wavenumber. A nonrobust embedded guided mode is characterized by an 
isolated point in the real frequency-wavenumber plane that lies on the complex dispersion relation. 

§5. Resonant scattering near a guided-mode frequency. This is the most important and 
interesting section of the paper. The complex perturbation analysis of [10j is applied to transmission 
anomalies for the discrete model. We extend it to capture the singular behavior of the transmitted 



energy at a tangent bifurcation of guided modes (Theorem 18), in which the bifurcation parameter 
is a constant of coupling between the ambient lattice and the waveguide. We also analyze the 
accompanying resonant amplification. 

2 The Discrete Model 

We have chosen to analyze a model in which the ambient space and the wave guide can first be 
described as separate systems in their own right, which are then coupled together through simple 
coupling constants along a line. 

2.1 The Ambient Planar Lattice 

The ambient space is a planar (two-dimensional) lattice of beads of mass 1 located at the integer 
points Z 2 in M 2 that are connected by springs of strength 1. The internal dynamics are given by a 
Schrodinger-type equation 

V = -Xhv, (2-1) 
where y = {y m n} £ ^ 2 (Z 2 ) =: H 2 with (m,n) 6 Z 2 and — Q 2 is the discrete uniform Laplacian 

(^2y)mn = -(y(m-l)n + V{m+l)n + Vm(n-1) + Vm(n+1) ~ ^Vmnj- (2-2) 

The spatial part u = {u mn } of a harmonic solution y mn (t) = e~ Mt u mn satisfies 

(Q 2 - oj)u mn = 0. (2.3) 

The solutions of this equation are generalized eigenfunctions of O2, and the simplest of these are 
the plane waves 

e 2«(me+n0) ) w = 4-2cos(27r0)-2cos(27n/>). (2.4) 

This relation between u and (9, (f>) is the dispersion relation for the free 2D lattice. 

Through the (inverse) Fourier transform, each element of H2 is expressed as an integral super- 
position of these eigenfunctions 
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The operator F is unitary. The bounded operator can be written in terms of the shift operators 
on H 2 , 

(S r u) rnn = U( m -i) n , (S U) mn = U( m _|-x)nj 
(S u u) rnn = U m (n—l^, (S u) mn = U m (n+l)i 

fl 2 = -(S l + S r + S u + S d -4I). (2.5) 
The shift operators become multiplication operators under the Fourier transform, 



TS*F- l f{e,<l>) = e- 2 " e /(M), 

FS l F- 1 f{e,<i>) = e 2 ™ e f{0,4>), 
FS n F- l f{9, (/)) = e~ 2 ^f(e, <j>), 
FS d F- l f{e,<f)) = e 2 ™ 4 



(2.6) 



and we thus obtain a spectral representation of O2, in which the value of the multiplication operator 
at (9, (f>) is the frequency given by the dispersion relation for plane waves above (2.4). 

(Fn 2 F- 1 )f(9, 0) = (4 - 2 cos(27r#) - 2 cos(2tt0))/(0, cp). (2.7) 

The spectrum of O2 is the range of this multiplier, [0, 8]. 

2.2 The Periodic Waveguide 

Our periodic waveguide is an infinite sequence of beads connected by springs. The internal dynamics 
in the Hilbert space H\ = £ 2 (Z) are given by the equation 

Mi = -iAx, (2.8) 

where M is the positive mass operator defined by 

(Mx)j := M jXj , Mj > 0, (2.9) 

the internal operator A is minus the discrete nonuniform Laplacian 

(Ax)j := —kjXj + i + (kj + kj-%)xj — kj-iXj-i, (2-10) 

and both M and A are taken to be iV-periodic: 

M j+N = Mj and k j+N = kj for all j G Z. (2.11) 

By redefining the variable by the substitution x 1— > M _ 5x and denoting the operator M _ 5.AM _ 5 
by Oi, we reduce the equation to a simpler form: 

x = -iSlix. (2.12) 

Since A and M are self-adjoint, so is Ox, and it is represented by a tridiagonal matrix with periodic 
entries, 

(n lX )j = - . kj x j+1 + ^X^ Xj - . kj ~ l Xj-i. (2.13) 
3 y/MjM j+l 3 Mj 3 yJMjMj-i 3 
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Since commutes with the shift operator S, 

(Sx)j = x j+N , 



(2.14) 



we can obtain by the Floquet theory the generalized eigenfunctions of by examining those of S. 
Since 5 is unitary, its generalized eigenfunctions {xj} are characterized by the pseudo-periodic 
condition 



Xj+N 



"2 < K — 2- 



(2.15) 



Let us denote by V K the ^-dimensional space of solutions, which is spanned by the vectors (i = 
1,...,N): 

pW = (...,e^^,0,...,0,^,0,...,0,e^,0,...,0,e^,...). (2.16) 

-N+i l N+e 2N+£ 



(k) 

With respect to this basis, the restriction of f2i to V K is represented by the Floquet matrix 

/ (fci+fcjy) -fci q ... q -fcjy c 2mK\^ 
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(2.17) 



Incident wave {a^} 




Figure 1: A schematic diagram of the discrete model. 



2.3 The Coupled System 

Let us couple the systems (Hi, Oi) and (H2, ^2) in a simple way by introducing a periodic sequence 
of constants 7 n with ^ n +N = 7n that couple x n to u$ n . This is achieved by the coupling operator 
r : H2 — > Hi defined through 

T(E 0n ) = 7„e n , T(E mn ) = 0, if m ^ 0, 
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in which {e n } ne z and {£ mn } m ,nez are the standard orthonormal Hilbert-space bases for Hi and 
H2, respectively. The adjoint : Hi — > H2 of T is 

^(e n ) = %E 0n . 

The internal dynamics in Hi and H2, together with the coupling between them, define a lossless 
oscillatory dynamical system (T~L,£l) in the Hilbert space 



U = Hi® H 2 , 

where f2 has the following form with respect to this decomposition 



n 



ni r 
rt ^ 2 



(2.18) 



(2.19) 



and the dynamics in % are given by 



The assumption of a harmonic field with circular frequency u>, x = e lujt z and y = e lujt u, leads to 
the eigenvalue problem 
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= -in 


X 
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y. 



n 



which is equivalent to the coupled system 



z 




z 




= UJ 




u 




u 



UlZ n = (Qiz) n + {Tu) n , 

uu mn = (r^z) mn + (^2^) 



(2.20) 



(2.21) 
(2.22) 



Because of the periodicity of the waveguide, the operator commutes with translation by N 
lattice points in the n variable, that is, u mn i-> u mtn -N- By the Floquet-Bloch theory, O is a direct 
integral of pseudo-periodic operators f2 K , 



n 



1 -} 

2'2> 



r2 K dnc, 



which are defined by the restriction of to the functions u mn that satisfy the pseudo-periodic 
condition 

lAm^+N — 6 U mn . 

2.4 Dynamics projected onto the waveguide 

If the conservative system (7i, £1) is projected onto Hi, the result is the dissipative system 

/>oo 

z = -m x z - \ Te- irn2 T ] z(t - r)dT , z(t) £ Hi. (2.23) 
Jo 

One can ask the question: how much of the original system (l~L,£l) can be reconstructed from the 



dynamics projected onto Hi, that is, from equation (2.23) alone; more viscerally: which motions 



of the ambient lattice can be detected by an observer living in the waveguide? Equivalently, one 
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could ask which motions originating in the ambient lattice can disturb the waveguide. Figotin and 
Schenker [5 J prove that a dissipative system of the form (2.23) admits a conservative extension 
(7i,Q), with dynamics v = —iVtv (v(t) £ H), that is unique up to Hilbert-space isomorphism. The 
unique extension of the projection of (T~L, Q) onto H\ can be realized as a unique subsystem of the 
original system f2). The construction of this subsystem is given by Theorem 9 in [6], 



U 

n = 



-H X 

J2i€ 



H 2 , H2 = On 2 (^(Hi)), 



n 



2\h 2 , 



in which the subspace H 2 of H 2 is the orbit of the image rt(#i) in H2 under the action of tt 2 . The 
orbit is defined by 

Definition 1 (orbit). Let Q be a self-adjoint operator in a Hilbert space H and S a subset of vectors 
in H. Then we define the closed orbit (or simply orbit) Oq(S) of S under action of VL by 



Oq(S) = closure of span {ip(Q)w : tp G C C (R), w E S}, 



(2.24) 



where C C (M) is the space of continuous complex-valued functions on M. with compact support. If H' 
is a subspace of H such that Oq(H') = H' , then H 1 is said to be invariant with respect to O or 
simply Q-invariant. 




Figure 2: Level sets of the multiplication operator J r Q 2 J r 1 : f(9, 4>) ^ (4 — 2 cos(27T0) — 2 cos(2ncj)))f(d, <j>). 

The following theorem says that the component of H 2 in the system (■%, Q) that is reconstructible 
from the dynamics projected onto Hi is the space of motions that are symmetric with respect to 
the variable m, in other words, motions that are anti-symmetric with respect to the line of coupling 
of the waveguide cannot excite the waveguide. 

Theorem 2 (part of H2 determined by Hi). The orbit ofT'(H\) in H 2 under the action of Q 2 is 

On 2 (F^(Hi)) = {{u mn } £ H 2 : u mn = u- mn \/m,n £ Z} . 

Proof. As we have seen, is represented on L 2 ([— l / 2 , V2] 2 ) through the Fourier transform J- by 
the multiplication operator T g f = gf, where 

g (e, <j))=A-2 cos(2vr(9) - 2 cos(2vr0). 

Furthermore, T maps the image 

V\Hi) = {{u mn } £H 2 :u mn = for m / 0} 
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to the subspace S of L 2 {[— / 2 , V2] 2 ) of functions that depend only on <j), 

S = H^(Hx)) = {he L 2 ([- % y 2 ] 2 ) : h(0, <t>) = hxid) for some h x E L 2 ([- % %})}. 

Since the space of functions u mn in i/ 2 that are symmetric in m is mapped by T onto the space of 
functions in L 2 ([— y 2 , /a] 2 ) that are symmetric in 0, it suffices to prove that 

closure of span ty(g(9, 0)M0) E L 2 ([-y 2) V 2 ] 2 ) : ^ E C(R),/ti E L 2 ([-V 2 , y 2 ])} 

= {/ e l 2 ([- y 2 , y 2 ] 2 ) : ho, <t>) = f(-e, $ vm e [- y 2 , y 2 ]}. 

Because 5 is symmetric in 9, all the functions in the space on the left-hand side of this equality are 
also symmetric in 9. Thus, the theorem will be proved by showing that 

closure of span {i>(g(9, <t>))hi{4>) E L 2 ([0, l / 2 ] x [- % %]) : $ E C(R), h x E L 2 ([- y 2 , y 2 ])} 

= l 2 ([o, y 2 ]x[-y 2 , y 2 ]). (2.25) 

Define 

s = : ^ e c(k),^ e c([-y 2 , y 2 ])} c c([o, y 2 ] x [-%, y 2 ]), 

21 = span(E). 

The set 21 is a complex algebra that is closed under conjugation, and it is therefore the algebra 
generated by the set E. By taking ifi and hi to be constant in the definition of E, we find that E 
contains the constant functions. To see that E separates points, let (6*i, <f>i) and (0 2 , </> 2 ) be distinct 
points in [0, x / 2 ] x [— y 2 , y 2 ]. If ^ 7^ 2 , then the function (0, 0) i-)- in E, obtained by taking tp 
to be unity and h\{4>) = <p separates these points. If 9\ ^9%, then, from the definition of g, we see 
that the function g(9,(f>) in E, obtained by setting ^(7) = 7 and h\{4>) = 1 separates the points 
(see Fig. [2]). By the Stone- Weierstrafi Theorem, 21 is dense in C([0, 72] x [— y 2 , y 2 ]) in the uniform 
norm. 

Since C([0, y 2 ] x [— x / 2 , x / 2 ]) is a dense subset of L 2 ([0, x / 2 ] x [— y 2 , y 2 ]) in the L 2 norm and the 
L 2 norm is bounded by the uniform norm, we conclude that 21 is dense in L 2 ([0, x / 2 ] x [— y 2 , y 2 ]). 



But, since 21 is a subset of the set on the left-hand side of (2.25), we have proved (2.25) and thus 
established the theorem. I 

3 The time- harmonic scattering problem 

We shall assume the harmonic time-dependent factor exp(— iujt) from now on. 
3.1 Spatial Fourier Harmonics 



Let us consider pseudo-periodic solutions to the problem (2.21), (2.22) with Bloch wave number k 



in the n-direction, which is the direction of the line of coupling. This means that 

Z n+N — E Z n , 
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Such solutions have finite Fourier representations: 

N-l 
£=0 



N 



u,, 



u,. 



N-l 

1=0 
N-l 



+ aJe- 27Tieem )e 2ni ^ n , m < 0, 



+ bp e 



-2iri9im\ 2iri(pgn 



m > 0, 



in which fa = (k + £)/N and 9i = 9i(k,oj) is the m-component of the wavevector determined by 
the dispersion relation (2.4) for the operator f?2; 

UJ = 4 - 2 C0S(27T^) - 2 COS(27T0£). 

Those values of £ for which 9% is real correspond to propagating Fourier harmonics (diffractive or- 
ders), and those values of £ for which 9e is imaginary correspond to exponential harmonics (evanes- 
cent orders). In the former case, we take 9g > 0, and in the latter, we take iOi < 0. These cases are 
separated by the case of a linear harmonic 9i = 0. 

Because of the periodicity of the structure, each pseudo-periodic function u mn is characterized 
by a minimal Bloch wave vector k lying in the first Brillouin zone k S [— 72, 72 )■ The region 
[— 72, 72] x [0,8] in (ac, w)-space is divided into sub-regions according to the number of propagating 
Fourier harmonics. For a given pair (k, oS), let ^3 be the set of propagating harmonics, 

<P = <p(/e, u) = {£ : Re(^) ^ 0}. (3.26) 

Figures [3] and |3] show the (k,uj) regions defined by the order of this set as a function of n and oj. 
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Figure 3: The diagram of |^3| for N — 2 (left) and N = 3 (right). The integers 0, 1, 2, and 3 represent the number 
of propagating harmonics. 

In the problem of scattering of traveling waves incident upon the waveguide from left and 
right, we must exclude exponential or linear growth of {u mn } in the ambient lattice as \m\ — > 00. 
Moreover, the energy of the scattered, or diffracted, field must be directed away from the scatter, 
that is, it must be outgoing. This notion is made precise by the following definition. 
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Figure 4: The diagram of |*p| for N = 9 (left) and N = 10 (right). 

Definition 3. (outgoing and incoming) A complex-valued function {u mn } is said to be outgoing if 
there are numbers {o^}* and {bi\^~Q such that 

N-l 

u mn = ^ a^-^ ie " m e 2 ^ n , m < 0, (3.27) 
e=o 

N-l 

Umn = b/,e 2 ™ e * m e 2 ^ n , m > 0. (3.28) 
The function {u mn } is said to be incoming if it admits the expansions 

N-l 

Umn = a e e 2nieem e 27Ti ^ n , m < 0, (3.29) 

N-l 

Umn = J2 b e e- 27Tieem e 27Ti4>en , m > 0. (3.30) 

1=0 

The form of the total field is the sum of an incident field and a scattered one and has the form 

N-l 

Umn = af c e 2nieem e 2ni ^ n + ^ ^-s^m^^n m < ^ (3 31 ) 

N-l 

in which the left-travelling second term for m < consists of the sum of the scattered field to the 
left of the waveguide and the incident field from the right with coefficients 6™°; the field for m > 
is understood analogously. 

Problem 4. (Scattering problem, P sc ) Given the coefficients {a™ c } and {^ nc } of an incident field, 
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find a pair of functions (z = {z n }, u = {u mn }) that satisfies the following conditions: 

UJZ n = (Qiz) n + (Tu) n , (3.33) 

WMmn = (^Z) mn + {VL 2 u) m n, (3.34) 

(z, u) are n-pseudoperiodic in n, (3.35) 

u = u mc + u sc , with u sc outgoing, (3.36) 

in which 

uZ n = [af c e 2widem + bf c e - 27Ti8em e 27ri ^ n . (3.37) 
3.2 Conservation of Energy 

Since we seek K-pseudo-periodic fields, the scattering problem can restricted to a strip 1Z containing 
one period with boundary in the variable n, 

K = {{m, n) eZ 2 : -oo < m < oo, < n < N}. (3.38) 

The coupled system admits a law of conservation of energy: The total time-harmonic flux ( 3.42| ) 
out of a truncated region [mi, 777.2] x [1,N] of the strip 1Z vanishes; this is stated in the following 
theorem. 

Theorem 5. Let the frequency u> and wavenumber k be real. If the pair (z n ,u mn ) satisfies the 
coupled system (2.21\2.22 ) and has the form 



N-l 



u,, 



u,, 



£ Qe 2 "*", 777 = 0, (3.39) 
=0 

N-l 

^2 (aje~ 2meem + aje 2meem )e 27Ti ^ n , m < 0, (3.40) 
e=o 

N-l 

Y (bje~ 27rieem + b+ e 2 ^ m )e 27Ti ^ n , m > 0, (3.41) 



then 



Im \ y2(uu x ) mir ^j - Im ^^(uu x ) m2 n \ = (3.42) 

Yl i(\ b t\ 2 + K~| 2 ) - (\4\ 2 + I&7I 2 )] sin(27r^) = 0, (3.43) 

where u x is the forward difference of u in the variable m (see the Appendix). The fluxes on the 
upper and lower boundaries cancel identically by pseudo-periodicity. 

Proof: We need the following summation-by-parts formula (see the Appendix): 

^/-r, \ k N z N Z N+ i z N k z , Z\ z s 

}2{z^iz) n = -==(—== -==) + —=(- 



n =i ^fM N ' y JM^ 1 VMn VMo y/Mi VW 

N _ - \ / X ^ 3 - 44 ) 

V - ^ , / z n Z n —l \ I Z n Z n —l 

n=l 
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Upon multiplying (2.21) by z n and summing over one period, the condition of pseudo-periodicity 
for z n and the above formula yield 



N 



n=l 



Zn-1 



N 



(3.45) 



n=l 



Similarly, multiplying (2.22) by u mn and using (6.119) we obtain 

N 



N m 2 



N 



N 



° = X, (w|n mn | 2 -|V_u mn | 2 ) + ^(™ x ) m2n -^(nu x ) min -^ 



UQnJnZn- 



(3.46) 



n=lm=mi rt=l n=l n=l 

The boundary values at n = and n = N cancel because of the pseudo-periodicity of u Ti 



Adding (3.45) with (3.46) and taking the imaginary part of the sum leads to the condition (3.42). 
Calculation of (3.43) is straightforward. I 



3.3 Formulation in Terms of Fourier Coefficients 

The scattering problem can be reduced to a system of equations for the Fourier coefficients: 

'jv-1 , -„ N-i , ■„ 

N-l „ JV-1 , „ 

< £ (aJe 2 ^-b+e- 2 ^-%c i )e^ ln = £ (bje 2 ™ 6 * -a+e" 2 ^) e^r n , (3.47) 
or in the matrix form: 

BA^ = ^, (3.48) 

where B is a 3Ax3A matrix, the vector F contains the coefficients {a/} and {67} of the source field, 

and the vector represents the coefficients {aj}, {&/}, of the outgoing field and the coefficients 
{q} of the field in the waveguide. 



3.4 Solution of the Scattering Problem 

To prove that the scattering problem always has a solution, it is convenient to work with its 
variational form. For this purpose we introduce artificial boundaries at m = —Ai and m = M. At 
these boundaries, the outgoing condition is enforced through an associated Dirichlet-to-Neumann 
operator T, which acts on traces on the boundaries m = of functions in the pseudo-periodic 
space 7i K (R,), 

U K (K) = {{z, u) e H(K) : z N = e 2niK z , u mN = e 2 ™u m0 }, (3.49) 

and is defined through the finite Fourier transform as follows. For any function v = {v n } n = , let 
Vg be the I th Fourier coefficient of the function {v n e~ 2mKn ^ N }^~Q , that is, 

JV-l 
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Then the map T is defined through 



(Tt/)? = (1 



(3.50) 



The operator T characterizes the normal forward differences of an outgoing function on the bound- 
ary m = =fA^ of the truncated domain in terms of its values there, 



(d v u + Tu)±Mn = for u outgoing, 



(3.51) 



where 



{d v u)-Mn = -U~Mn + U(-M-l)n, (3.52) 
{d v u)Mn = u (M+l)n ~ UMn- (3.53) 

Then using the decomposition u = u sc + -u mc of the solution to the scattering problem P sc we obtain 



d v u + Tu 



d u u mc 



+ Tu mc 

2 E C 1 - cos(2vr^))a+e 27r ^(-- M )e 27r ^ n , m 



-M, 



(3.54) 



2 E (1 - cos(27r0*))&7e- 



-2-Kido M 2irid>?n 



m 



M. 



Thus we are led to the following problem set in the bounded domain TZ M of Z 2 : 

n M = [-M,M] x [0,Nj. 

Problem 6. (Scattering problem reduced to a bounded domain, Pjfi) Find (z,u) in %{1Z 
that 



ujz n = (p,xz) n + (Tu) n for < n < N, 

uu mn = (T^z) mn + (0, 2 u) mn for (m,n) £ TZ M , 

(z, u) are K-pseudoperiodic in n, 



3 v u + Tu = d u u mc + Tu mc on m 



(3.55) 
M ) such 

(3.56) 
(3.57) 
(3.58) 
(3.59) 



Problems P sc and are equivalent in the sense of the following theorem. 



Theorem 7. If (z,u) is a solution of P sc such that (z,u) = (z,u)\ r m £ 7-L(R M ), then (z,u) is a 
solution of Pj^. Conversely, if(z,u) is a solution of Pj^, it can be extended uniquely to a solution 
(z,u) ofP sc . 



3.59 



Proof: The first part of the theorem holds because the condition (3.59) is equivalent to (3.36) 
Conversely, if (z,u) is a solution of Pffi, then by 



the difference n s 



u 



satisfies the 



Dirichlet-Neumann relation (3.51) that characterizes outgoing fields and can therefore be extended 
to an outgoing field u sc . The field u = u mc + u sc , by the definition of its parts, satisfies ( 3.33|3.34 ) 
outside of K M . I 

A variational form of the scattering problem, which is analogous to the weak formulation for 
partial differential equations, is obtained from the summation by parts formulas in the Appendix. 
With the notation d y and dy for the forward and backward differences in the n-variable, one can 

write Qi = — Ma dyKd y M?- , where (Kv) n = k n v n , and use the substitution v h-> Kd y M^z and 
. i 



w i y Mz-w in (6.116) to obtain the first equation below. The second is obtained by using equation 
(6.119). We use the notation V_ = (dx,dy). 

It is straightforward to prove that problems Pj^ and PJ^j. are equivalent. 
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Problem 8. (Scattering Problem, variational form, PJ^-J Find a function (z,u) G rl K (R M ) such 
that 

N-i i _ i _ TV 

]T {Kd y M 2 Z ) n (d y M 2w) n + ((Tu) n w n - (ujz) n w n ) = 0, 

n=0 n=l 

A 1 " M N _ 

E J] (wm)-V_-u V_n-(rt 2 ;)w) mn )-^(v(-A4-i)n(Tu)_ 7 v!n-(wT'ii)^(„) 

n=lm=—A4 n=l 

N N 

E (^(->l-i)n(^n inc +rn inc ))_^ n - £ (^u inc +TV nc ))A< n . 



(3.60) 



n=l 



n=l 



/or a// u) G 'H K (1Z- M ) and n-pseudoperiodic v = {v n }. 



Theorem 9. (Equivalence of P^ and PJ^.) If(z,u) 6 % K (7£ -A/! ) satisfies the scattering problem Pj^, 
then (z,u) satisfies PJ^.. Conversely, if (z,u) satisfies PJ^. for any (w,v) G % K ,(1Z M ), then (z,u) 
satisfies Pj^ also. 

The scattering problem always has a solution, even if it is not unique. Non-uniqueness occurs 
when the structure supports a guided mode, as we discuss in the next section. The reason for the 
existence of a solution of the scattering problem in the presence of guided modes lies in the orthog- 
onality of guided modes to incident plane waves: the former possess only evanescent harmonics, 
while the latter possess only propagating harmonics. This fact will be important in the analysis 
of resonant amplitude enhancement, and it has its analog in continuous problems of scattering of 
waves by open periodic waveguides [H Thm. 3.1] 



Theorem 10. The problem (3.60) always has a solution. 



Proof: Let us rewrite (3.60) in the concise inner product form 

(AY, V) = (F, V) for all V, 



(3.61) 



where Y = (z,u), F = (0, /) and V = (w, v). We use the Fredholm alternative, namely, that ( 3.61| ) 
has a solution (z, u) if and only if (F, V) = for all V G Null(^4*), or, in other words, 



(F, V) = for all V such that (AY, V) = for all Y. 



(3.62) 



Any function V satisfying (AY, V) = for all Y satisfies (AV, V) = as well. By Theorem [5] it 
follows that v contains only evanescent harmonics (or linear ones for threshold values of (k,cj)) in 
its Fourier series, that is 



— — 2iri0/m „2md>m 

E,.+ „2wi9em Jlm&in 
_ u e e e 



for m < 
for m > 



(3.63) 
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where = 1, . . . , iV \ Using the orthogonality of the Fourier harmonics, we obtain 
(F,V) = (f,v) 

N 

= ~ 2 E ((E vJe 2vie ^- M -^e- 2 ^ n ){Y, (1 - cos(27r^/))a+e 2 ^'(-^)e 27r ^ m )) 

n=l t'e<$ 
N 

~ 2 E ((E vj-e- 2m6eM e- 27ri ^ n )( EC 1 - cos(27r^))6^e- 27ri ^'- M e 27r ^ m )) 
= - 2 E_ ((^e 2 ^(-^- 1 )(l-cos(27r^))a+e 2 ^'(-^))(^e 27ri (^'"^) n )) 

<g<p n=l 

AT 

-2E_ ((^e- 2 ^ e ^(l-cos(27r^))6^e- 2 ^'^)(£e 2 ^(^'-^) n )) = 0. 

<g<p n=l 



Therefore there exists a solution (.z,u) to the problem (3.60). I 



4 Guided modes 

A guided mode is a nontrivial solution of problem P sc in which the incident field is set to zero. 

Problem 11. (Guided mode problem, P9 m ) Find a pair of functions {z = {z n },u = {u m n}) that 
satisfies the following conditions: 

U)Z n = (ttiz) n + (Tu) n , (4.65) 

= (T^z)mn + (fi2«)mn, (4.66) 

(z, u) are n-pseudoperiodic in n, (4.67) 

u is outgoing. (4.68) 



Because of conservation of energy relation (3.43), a generalized guided mode supported by the 



structure at a real pair (k,cj) possesses no propagating harmonics and is therefore exponentially 
decaying as \m\ — > oo. It may be called a true guided mode. 

In a real (k,uj) region in which all harmonics are evanescent (see the region labelled "0" in 
Fig. [3]), there are, for appropriate values of Mi and 7^, dispersion curves defining the locus of 
(k, w)-pairs that support a guided mode. These are robust in the sense that, as k is perturbed, the 
guided mode persists albeit at a different frequency. Physically speaking, energy cannot radiate 
away from the waveguide because there are no propagating harmonics available for transporting 
the energy (the energy flux of all evanescent harmonics is along the waveguide). 

The situation is different in a (k, w)-region of, say, one propagating harmonic (the region labeled 
"1" in Fig. [3]), as it typically contains no dispersion curves for guided modes. Nevertheless, a 
structure may support a guided mode at an isolated (/i,w)-pair in this region. Such a guided 
mode is nonrobust with respect to perturbations of k. The physical idea is that the propagating 



harmonics, which appear with nonzero coefficients in a typical solution (3.39-3.41) are carriers of 
incoming and outgoing radiation. But special conditions that allow these coefficients to vanish at 
some frequency, thereby creating a guided mode, may be arranged by tuning the structure (masses 
and coupling constants) and the wavenumber to specific parameters. A perturbation of k from the 
value kq that supports the guided mode will destroy these special conditions. In physical terms, one 
often describes the destruction of a guided mode as the coupling or interaction of the mode with 
radiative fields. It is this interaction that causes resonant scattering behavior and transmission 
anomalies. 
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4.1 Generalized guided modes 

A rigorous analysis of resonance near nonrobust guided modes requires an extension of the scattering 
problem (P sc ) and the problem of guided modes (P gm ) to complex k and oj in a vicinity of the 
guided mode pair (kq,ljq). There arises a complex dispersion relation D(k,uj) = in C 2 describing 
the locus of generalized guided modes, given by the zero set of the determinant of the matrix B in 



(3.48). From this point of view, an isolated guided-mode pair in real (k, w)-space is the intersection 
in C between the real plane and the dispersion relation. We call solutions of Problem P gm for 
complex k or w "generalized guided modes" ; they are foundational to the theory of leaky modes, 
as discussed, for example, in [71 [T3"l IT4"1 121] . 

Let us consider real values of k and examine how the nature of the Fourier harmonics changes 
when u is allowed to assume a small imaginary part. Suppose that 0i > for some real u = ur, 
that is, the £ t]l harmonic is propagating. Then for to = ojr + itaj, where uj is small, 0g also attains 
a small imaginary part, 0g = Of' + i0l. The dispersion relation gives 

2cos(2vr((9f + i0{)) = A-oj R -iuj I - 2 cos(2vr0/;). 

Taking the imaginary part, one finds 

sm(2ir0?)(e 2 <-e- 2 * 6 l)=cj I . 



Thus, if loj is negative (and sufficiently small), 0\ is also negative. This means that the outgoing 
Fourier harmonic e 2m (<t>i n +Qi\ m \) e —iut Qi eca y S [ n time but grows in space as \m\ —> oo whereas the 
incoming one decays in space and time. Conversely, if uj is a small positive number, then 0\ > 
and the £ th outgoing harmonic grows in time and decays in space and the incoming one grows in 
space and time. The evanescent (resp. growing) harmonics, for which 0g < 0, remain evanescent 
(resp. growing) under small imaginary perturbations of oj. 

The following theorem is the discrete analog of Theorem 5.2 in [18] for photonic crystal slabs, 
which states that generalized guided modes occur only for Im(u;) < and that such a mode is a 
true evanescent one if and only if Im(w) = 0. 

Theorem 12. Suppose that ({z n },{u mn }) is a nontrivial solution to the homogeneous (sourceless) 
problem P9 m . Then Im(u;) < 0. In addition, \u mn \ — > as \m\ — > oo if and only i/Im(a;) = 0. 

Proof: Adding ( 3.45[ ) and (3.46) and taking the imaginary part yields 



TV 7712 N N 

lm(u)(^2 \ z n\ 2 + I"™™! 2 ) = Im(y^(mt x ) min - ^(uu x ) m2n ) . (4.69) 

77 = 1 71=1 777=7771 77=1 77 = 1 



If the field decays as \m\ —> oo, then the right-hand side of (4.69) tends to zero as —mi and 
tend to oo, and thus the left-hand side vanishes. Since the field is nontrivial, we obtain Im(w) = 0. 
Conversely, if Im(w) = 0, then Theorem [5] applies and we find that, since the coefficients a\ and 
bj {£ € ^3) of the incident field vanish, so also must and b\ for £ G <p. Thus u mn decays 
exponentially as \m\ —> 00. 

If Im(u) > 0, then, as we have mentioned, the generalized outgoing Fourier harmonics of the 
field decay as \m\ — > 00. Since there are no incoming harmonics by assumption, the right-hand side 



of (4.69) decays as m — > 00, so the left-hand side vanishes, and we must have Im(cj) = 0. We 



conclude that, for all solutions of the homogeneous problem P§ m , it is necessary that Im(w) < 0. I 
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4.2 The spectrum of fi K 

As we have mentioned earlier, the operator Q is a direct integral of the self-adjoint operators Q K . 
Since the latter acts on the space of K-pseudo-periodic functions u mn , its domain may be taken to be 
1 2 (1Z) with the pseudo-periodic boundary condition u m N = e 2mK u m o- The spectrum of f2 K consists 
of a continuous part and a set of eigenvalues. The continuous spectrum is the ^-coordinates of the 
intersection of the region of at least one propagating harmonic in Figs. [3] and [4] with the vertical 
line at fixed k. The eigenvalues are the real values of oj for which B has a nullspace, that is, the 
frequencies that support a true guided mode for the given value of k. These frequencies may be 
either in the region of no propagating harmonics or embedded in the continuous spectrum. 

In the context of spectral theory, an isolated point (kq, ujq) of the dispersion relation D(n, u) = 0, 
within a region of at least one propagating harmonic, corresponds to an embedded eigenvalue of the 
operator Q K that dissolves into the continuous spectrum as k is perturbed from kq. The associated 
destruction of the guided mode is associated with transmission resonance. This phenomenon is 
akin to the quantum- mechanical resonances of the noble gases, in which an embedded bound state 
of the idealized atom with no interaction between the electrons is destroyed when this interaction 
is initiated, @], [HI §XXII.6]. 



4.3 Existence of Guided Modes 

Let us consider the existence of guided modes in the case of period N = 2, which is a minimal model 
for the phenomenon of anomalous transmission. Depending on the values of k and w, there are 
either two, one, or no propagating harmonics, as depicted in Fig. [3j We are interested in regions I 
and II, in which there is exactly one propagating and one evanescent harmonic. It is this region 
in which one may encounter a nonrobust guided mode and associated scattering resonance. 

In the region /, the exponent 9q corresponds to a propagating harmonic, whereas 6\ corresponds 
to an evanescent one. The coefficients of the propagating harmonic in the Fourier decomposition 
of the solution to the guided mode problem (P gm ) are forced to be zero. The corresponding 
unknown vector in equation (3.48) is X = (<Zg = 0, b^ = 0, cp, aj, c\ ). The condition that the 
values cq, a7 , bt, and c\ not vanish simultaneously in the system (3.47) yields two equations that 



characterize a guided mode 



(7i ~ 7o) 
(7o + 7i) 



(fco + fci)( 



1 1 

Mi ~ Mo 
7o7i(7o + 7i) 



(70 + 7i)£sin(27r0i) 



2i sin(7TK) 

+ , (kg - h 

+ 2u + (k + k 1 )(- 



1 

Mo 



1 

Wi 



2 cos(7tk) t 
VMqMi ' 



(4.70) 



(71 ~ 7o) 
(7o + 7i) 



2w + (/c + /ci)( 



2 cos(-7tk) 



+ 



VM Mi 
7o7i(7i - 7o) 



1 

Mo 



1 

Mi 



(7o + 7i)*-sin(27r0i) +(fc ° + fcl)( M 1 



M - 



+ 



2% sin(-7rK)(A;i — ko 
VM Mi 



0, (4.71) 



where sin(27r^i) = -y/l — (2 — | + cos(7tk)) 2 . For any real pair (ko,^o) satisfying ( 4.70|4.71 ) with 
no £ [ — 72L we must also ensure that oj§ G (2 — 2 cos(7tko), 2 + 2 cos(7r«o))) that is, that (kq, ujq) 
lies within region /. 
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Region // can be treated similarly. There, the roles of 6q and 0\ are switched and one seeks 
nonzero solutions for X = {clq , b$, cq, = 0, bf = 0, c\). 

Figures [5] and [7] show the locus of the solutions of ( 4.70|4.71~ ) in the (k,uj) plane for different 
choices of the structural parameters of the system. The points of intersection come in ±k pairs 
when the k- value of the corresponding guided mode is nonzero (Fig. [5]). Fig. [7] shows the case of a 
guided mode at k = kq = 0, which, as the system parameters are perturbed, either bifurcates into 
two modes at ±kq 7^ or disappears altogether. We will examine this bifurcation in more detail in 
Section [5j 

For certain structures, one can prove the absence of non-robust guided modes. 

Theorem 13. For the coupled system of period two with Mq 7^ M\, ko = k\ = k 7^ ; and 

±70 = 71 = 7 7^ 0, there is no guided-mode pair (k,u) in the subregion of [— 72, 72] x [0,8] that 
admits at least one propagating harmonic. 

Proof: Suppose the real pair (k,oj) admits a solution {u mn , z n } of Problem p9 m . By Theorem [T2| 
\umn\ —t" as \m\ —7- 00, and thus the solution contains only evanescent harmonics. Given that (n, oj) 
admits at least one propagating harmonic, we must have |*}}| = 1 (see Fig. 3, left). In the region 
/ defined by k G [—72, 72], to £ [2 — 2cos(7tk),2 + 2cos(7tk)], with propagating and evanescent 
harmonics corresponding to #0 £ K and 9\ E iM., respectively, the solution has the form 



U"mn — °"i e t 



m < 0, 



Umn = ft+ e 27ri0im e 2^ 1 n j m > Q) 

Zn = c e 2ni ^ on + cie 2?r ^ in , 



(4.72) 



with a, 







corresponding to the propagating harmonics. Using the restrictions on ki and 



7i given in the Theorem, the system (3.47) with zero right-hand side can be reduced to derive 
Co = c\ = and = bf = 0. Region iY is handled similarly. I 
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Figure 5: Example for Mo = 2, Ml = 1, ko = fci = 1, 70 = 1, and 71 = 7. Left: The intersection of the two 
relations ( 4.70|4.71 l (in red and blue) is the locus of guided modes, in this case (ko = ±0.0616, ojo = 0.9792). A region 
of one propagating harmonic is outlined in black. Right: The unfilled circles outline the region of one propagating 
harmonic and indicate the real part of the dispersion relation within this region. The pair (ko,u;o) of a guided mode 
in this region is represented as a solid black dot, as is the real dispersion relation for guided modes in the region of 
no propagating harmonics. 



One can construct guided modes for the system with period ./V greater than 2. By choosing M n , 
k n , and 7 n to be symmetric (about some index n if N is odd and about a point between some n 
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0.974 0.976 0.978 0.980 0.982 0.984 



Figure 6: Left: The transmission coefficient for Mo = 2, Mi = 1, ko = k\ = 1, 70 = 1, 71 = 7. There is a guided 
mode at (ko,wo) ~ (0.0616,0.9792). Right: Refined graphs of T vs. ui for various values of k near no- Observe that 
i\ 7^ in the formula of Theorem [16] for this example, so that the center of the resonance varies linearly as a function 
of k — Ko, whereas the width varies quadraticalfy. 



and n + 1 if ./V is even), one can construct anti-symmetric guided modes at k = in a (k, to) region 
that admits only one propagating harmonic. The idea is that, for k = 0, the system decouples into 
symmetric and anti-symmetric parts, and the single propagating harmonic is constant in n and 
therefore necessarily even. One then seeks solutions to the anti-symmetric guided-mode problem, 
in which propagating harmonics are automatically absent. This idea underlies behind the existence 
of non- robust guided modes in open electromagnetic or acoustic waveguides [201 E] • 

In the case of period three, for example, we can make the following specific assertion. 

Theorem 14. For the coupled system of period N = 3, 

1. if ji = 1 and ki = 1 for i = 0,1, 2; and < M\ = M 2 = M < 3^21 and M / M > 0, there 
is a guided mode at k = and ui G (0, 3); 

& if 7o = 7i = 72 = 7 and if there is a guided mode at k = and u G (0,3) then cq = 0. 
Additionally if Mi = Mi = M > 0, then the mode is antisymmetric about a line passing through 
the zeroth bead, that is c\ = — ci, oJ~ = — a% , and bf = — 6^7 

Proof: To prove part (1), the symmetry of the structure about a line passing through the zeroth 
mass allows the construction of an antisymmetric guided mode at k = about the same line. The 
coefficients Cj for the field in the waveguide satisfy 

c = 0, ci = -c 2 , (4.73) 

while the coefficients for the field in the ambient lattice satisfy 

a =b+ = 0, = -aj, b+ = -b+. (4.74) 



Taking into account these restrictions on the coefficients, the system (3.47) with vanishing right- 
hand side becomes 

d = 2m5"sin(27r6>i), (4.75) 
^sin^TrfliXw-IO- 1 = 0, 
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-0.4 -0.2 0.2 0.4 0.1 0.2 „ 0.3 0.4 O.J 

K K 

Figure 7: Example for a non-robust guided mode at n = with Mo = 2, Mi = 1, k = ki = 1, 70 = 1.029633513, 
and 71 = 7. Left: The intersection point of the two relations gives the parameters (/to — 0,ujo ~ 1) of the guided 
mode. Right: Real part of the dispersion relation in the region of one propagating harmonic. 



where sin (2tt0i) = iy (^y^) 2 — 1- As long as M < 2\/21, the third of these equations determines 
oj € (0, 3/M). The first two equations complete the determination of the Fourier coefficients, 
providing a one-parameter of family of guided modes. 

To prove part (2), one sets Oq = b 
vanishing right-hand side. I 



and deduces the relations ( 4.73|4.74 ) from (3.47) with 



5 Resonant scattering 

The analysis of resonant scattering at wavenumbers and frequencies near those of nonrobust guided 
modes is based on the complex-analytic connection between guided modes (for D(k, uj) = 0) and 
scattering states (for D(k,u) 7^ 0). This is achieved by scaling the incident field by an eigen- 
value £(k,uj) of B whose zero set coincides with the dispersion relation near a real guided-mode 
pair (/co, wo). A complex perturbation analysis of the the eigenvalue, the complex transmission 
coefficient, and the complex reflection coefficient, all of which vanish at (/cq,u;q), yields asymptotic 
formulas for transmission anomalies. 

The analysis below follows that of [19] for anomalous scattering by periodic dielectric slabs. 
There, the Bloch wavenumber of the nonrobust guided mode vanishes kq = 0, so that the mode is a 
periodic standing wave. The consequence of this is that the transmission anomaly, to linear order in 
re— /Co, remains centered about the frequency wo of the guided mode. The existence of standing modes 
whose (minimal) period is equal to that of the structure can be proved |20l [1] , but there seems 
to be no proof hitherto in the literature of the existence of guided modes at nonzero re. (Guided 
modes for re =f whose period is a multiple of that of the slab can be constructed, but these are 
described by a real dispersion relation and are therefore robust under perturbations of re). 

It turns out that one can show that our discrete model with period N = 2 admits truly traveling 
modes (reoT^O) for appropriate choices of the masses and coupling constants. We will see that the 
nonvanishing of reo coincides with a shifting, or detuning, of the central frequency of the resonance 
as re is perturbed from kq- This detuning is linear in re — kq, while the width of the resonance 



increases as (re — kq) ■ In section 5.2, we show how standing modes and traveling modes are 
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0.977885910 0.977885915 0.977X85920 0.9778X5925 0.977885930 0.977885935 



Figure 8: Left: The transmission coefficient for M = 2, Mi = 1, 70 = 1.0296335f3, 71 = 7. There is a guided mode 
at (/to, Wo) = (0,0.9778859328...). Right: with refining resolution for k. Observe that, for this example, l\ = in the 
formula of Theorem |16[ so that both the width and the center of the resonance vary quadratically as a function of 




Figure 9: The transmission coefficient for N = 3 with 70 = 71 = 72 = 1, fco = fci = = 1, Mo = 1, Mi = M2 = 2, 
so = 0, ujo = 1.191465768 (left), with refining resolution for k (right). 

connected through a structural parameter. Keeping all parameters fixed except 70, which controls 
the coupling of the even-indexed sites of the waveguide to the planar lattice, we show that, if 70 lies 
above a certain critical value, the system admits no guided modes, that a standing mode (kq = 0) is 
initiated at the critical coupling, and that this mode bifurcates into two guided modes traveling in 
opposite directions (kq = ±k) as 70 passes below the critical value. The behavior of the transmission 
coefficient is complicated near the point of bifurcation, and we give an asymptotic formula for it 
that incorporates k, u), and 70. 

5.1 Asymptotic Analysis of Transmission Near a Guided-Mode Frequency 

Nontrivial solutions of the sourceless problem = occur at values of k and co where the matrix 
B has a zero eigenvalue I = £(k,uj) = 0. The relation £(k,oj) = 0, or co = oj(k) when solved for oj, 
is a branch of the complex dispersion relation for generalized guided modes. We analyze states 
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that correspond to a simple zero eigenvalue £ (that is, having multiplicity 1) occuring at a real 
pair (ko,cjo) that is in the region I in Fig. [3] with one propagating harmonic corresponding to 9q. 
By of the analyticity of £{k,uj) and under the generic assumption that d£/dui 7^ at (kq, ujq), the 
Weierstrafi Preparation Theorem provides the following local form for the dispersion relation: 



£(k,uj)=0 uj = uj + £i(k - «o) + £ 2 {k - k ) 2 + 0(\k - K | a ) 



where £\ is real, and Im^2 > because Imw for real k cannot be positive due to Theorem 12 



Since £ is of multiplicity 1 near (fto, ujq), there is an analytic change-of-basis matrix C such that 
B has the form 

n = CSC- 1 = C^ Q ^C- 1 (5.76) 

where the analytic matrix IB has dimension (3N — 1) x (3N — 1) and a bounded analytic inverse. 
Let (p mn (K,Lo) be a given analytic source field, as a plane wave incident upon the waveguide from 
the left, 

_ 2ni(e m+<j>on) 



The source vector F in (3.48) is then determined by (3.47) and this choice of source field, and F 
can be decomposed into its resonant and nonresonant parts 

p = aCe l + C{Q,F 2 ), (5.77) 

where the complex scalar a = a{n,uj) and the vector F2 G C 3Ar_1 , are analytic, and the vector 
e\ = (1, 0, . . . , 0) GC 37V . We now scale this source by a constant multiple of £, c£(k, u)(p mn (K, uj), so 
that it vanishes on the dispersion relation near (kq,coo) (c / is to be fixed later), and solve 

MX = c£F. (5.78) 

The solution is 

X = caCex + ciC{Q,B~ 1 F 2 ). (5.79) 

The analytic vector X corresponds to a solution field ip mn (K,uj) that connects scattering states 
with guided modes. If £(k,lo) = 0, ip mn is a generalized guided mode, otherwise it is a scattering 
state. 

For (k, uj) near (ko,ujo) for which £(k,uj) 7^ the solution -i/w in the ambient lattice satisfies 
the asymptotic relations 

ipmn ~ be 27Ti6om e 2ni(l}on m^oo. 



By the conservation of energy relation (3.43), for real (k, uj), we have \£\ = \a\ + \b\ , which 
implies that £, a, and b have a common root at (kq,loq). In the following analysis we use the 
notation uj = uj — uiq and H = k — kq. 

The Weierstrafi preparation theorem for analytic functions of two variables provides the follow- 
ing forms for £, a, and b : 

£{k,u) = e l P^ + £ l K + £ 2 n 2 + 0{\k\ z )][1 + 0{\k\ + |w|)], 
a(n,uj) = e ip2 [uj + r{k + r 2 K 2 + O(|K| 3 )][r + 0(\k\ + \uj\)], 
b(n, uj) = e^ 3 [uj + t X K + t 2 K 2 + 0(\k\ 3 )} [t + 0(\k\ + \u\)}, 
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in which ro, and to are positive real numbers and we have chosen c so that the corresponding 
coefficient for £ is unity. Using these expressions, we expand the relation |^| 2 = \a\ 2 + |6| 2 for real 
(k, uj) to obtain the following relations among the coefficients: 



1 = r 2 , + tl (uj 2 term), 

£\ = rl\ri\ 2 + tl\ti\ 2 (k 2 term), 

i\ = rQRe(ri) + igRe^i) (ojk term), 

Re(£ 2 ) = rgRe(r 2 ) + 4Re(t 2 ) (ujk 2 term) 

hRe(£ 2 ) = rgRe(r 2 fi) + ^Re(t 2 ti) (k 3 term), 

|^| 2 +2£iRe(£3) = ^[|r 2 | 2 +2Re(rif3)]+^[|t2| 2 + 2Re(tit3)] (k 4 term). 



(5.80) 



Because of Theorem |12| which says that Imoo < if k is real and £(k, uj) = 0, we find that £\ must be 
real-valued and that Im£ 2 > 0. Because of the equations r 2 , + t\ = 1 and l\ = rQRe(ri) + ipRe^i), 
£i lies between Re(ri) and Re(ii). 



Theorem 15. £ 1 =ti = r 1 e R. 

Proof: Suppose r\ = r\ R + ir\j and t\ = ti R + itu, then it follows that 

h=r 2 r 1R + t 2 t 1R and i\ = r 2 r 2 R + t 2 t 2 R + r 2 r 2 u + t 2 t 2 u . 
By convexity, the first equality implies 

l\ < rlrf R + ifo\ R . 

This is consistent with the second equality if and only if r\i = tu = and r\ R = t\ R , which yields 

h = n = ti. I 

We show now how to obtain a formula that approximates the transmission anomalies. According 
to the above theorem we use the expansions for a and b including terms of the second order in 
that is 

£ = e ipi(Q + £{k + £ 2 k 2 + . . .)(1 + ci£5 + c 2 K + . . .), 

a = r e ip2 (u + £{k + t 2 n 2 + ...)(!+ piuj + p 2 K + .. .), (5.81) 
b = toe^ 3 (£ + £{k + r 2 K 2 + . . .)(1 + q x Zj + q 2 K + . . .). 

In the first factors, the higher-order terms are 0(|/5| 3 ), in the second, they are 0(k 2 + uj 2 ). The 
transmission cefficient T depends on the absolute value of the ratio b/a, 



\b\ \b\ \b/a\ 



\£\ yff+f v/i + lW 

and b/a has form 

b i *o (2 + t X K + t 2 K 2 + 0(\k\ 3 )) 



(5.82) 



a ro (uj + r{k + r 2 K 2 + 0(|/5| 3 )) 
in which p = p% — p 2 , 771 = q\ — p\ , r\ 2 = q 2 — p 2 . The approximation 

b to\uj + tlK + t 2 K 2 ' 



(l + J7ia5 + 772/c + 0(|K| + |o;|)), (5.83) 



a 



±\l + rju\ (5.84) 
ro\u) + riK + r 2 K z | 
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Figure 10: The transmission coefficient for Mo = 2, Mi = 1, 70 = 1, 71 = 7. Upper: The exact formula for 
the transmission coefficient. Lower: Approximation with second order term in re, rj w 0.767728, to ~ 0.3142988, 
ro « 0.94932. A guided mode is supported at the parameters («o,wo) ~ (0.0616,0.9792). 



yields the following approximation for the transmission coefficient 

2 tl\bJ + t 1 K + t 2 K 2 | 2 |l + T]Uj\ 2 



uj + r\K + \ + ^ol^ + + | |1 + 



(5.85) 



which agrees well with the exact formula (see Fig. 10). One can see on those graphs that a sharp 
resonance emanates from the guided-mode frequency ujq as the wave number k is perturbed from kq. 
The anomaly widens quadratically as a function of H and it is detuned linearly away from the guided 
mode frequency ujq, which indicates that t\ 7^ 0. This is formula generalizes that of [19j, where it 
was assumed that l\ = because the structure was symmetric with respect to a line perpendicular 
to the waveguide and the guided mode was a standing wave (kq = 0). 



One can prove that one commits an error of 0(\k\ + uj ) in the approximation (5.85). 



Theorem 16. If £(k,uj) has a root at (kq,ujq) G M 2 ; the partial derivatives of I, a, and b with 
respect to uj do not vanish at (kq,ujq); an d Im^2 ^ in the form (5.81), then the error in the 
approximation (5.85) is of order 0(\7t\ + uj 2 ) and the following approximation holds: 



T(k, uj) = 

as (k,uj) — > (0,0) in M 2 , where £i 



to 



\uj + i\K + t{k l \ 

\uj + £lK + e2K 2 \ 



\i + CM + o(\k\ + uj 2 



(5.86) 



qi - ci . 
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Proof: We shall prove only the formula (5.86); the error in (5.85) can be proved similarly. 

(£5 + l{k + t 2 H 2 + t 3 K 3 + ■■■) 



&(«, to) 



£(n, 



to 



(to + £ lK + £ 2 k 2 + £ 3 k 3 + 
Since £\ is real- valued and \m.£ 2 7^ 0, the denominator can be written as 



Denote e 



u+£ 1 K+£ 2 K 2 +£ 3 K 3 +0(\^) = (u+£ 1 K+£ 2 k 2 ) (l + l I K 
£ 3 k 3 + 0(\k\ 4 ) 



(5.87) 



(5.88) 



LO + £ 1 K + £ 2 K 2 
(g + £{K + t 2 K 2 + t 3 K 3 + ■■■] 

(00 + £{ii + ^£ 2 + £ 3 H 3 + ■■■ 
(to + £{k + t 2 K 2 ^ 



0(|k|). Using |e| < 1 we obtain 



(to + + t 2 K 2 + t 3 K 3 + ■ ■ ■ 



(to + £ik + £ 2 n 2 ) 



CJ + + £ 2 K 



(l-e + e 2 



(5.89) 



(w + ^lAt + £ 2 K 2 ; 



(l-e + - 



Thus the expression in absolute values on the right of (5.87) is 

;2\ 



(w + ^k + ^k 2 ) 

(E~3^) 



+ 



[l -€ + •••)(! + C1W + C2K + -") 

-(1 -e + •••)(! + Ciw + C2« + " 



(5.90) 



(w + £ik + £ 2 k 2 ) 

Again, because Im^ 7^ 0, the second term is 0(|re|), and we obtain the result. 
5.2 Bifurcation of Guided Modes and Resonance 

In this section, will see how the strength of the coupling between the waveguide and the ambient 
lattice acts as a tangent bifurcation parameter for the creation and splitting of a guided mode. We 
will study the simplest case of period 2 with ko = k\ = k, in which we fix all parameters except 
one of the coupling constants. When this constant is lowered to a specific value, a guided mode is 
created at re = kq = (Fig. [5]) and is thus a standing wave exponentially confined to the waveguide. 
When the constant is lowered further, the mode splits into two guided modes at re = ±kq 7^ 
(Fig. [7]) traveling in opposite directions along the waveguide. 

Such a tangent bifurcation connects the case of £\ = in the transmission coefficient in Theorem 
16 to the case of £\ 7^ 0. Indeed, for a standing mode (reo = 0), re = k in the theorem and the 
symmetry of the transmission coefficient with respect to re implies that £\ = 0. In this case, the 
transmission formula shows that both the center of the anomaly and the distance between the peak 
and dip vary quadratically in re. On the other hand, when reo 7^ and the mode is traveling, we 
typically have 4/0 and thus there is a detuning of the central frequency of the anomaly away 
from the frequency too of the guided mode; the detuning is related linearly to re (to leading order), 
while the anomaly widens still only at a quadratic rate. 

We will need the following technical lemma. 

Lemma 17. Suppose that, for N = 2 and for fixed real values of Mo, M\, ko = k\ = k, 70 and 71, 
there is a unique real pair (kq,coo) i n an open set U of the real (re, to) -region I (Fig. [3]) of one 
propagating harmonic that admits a true guided mode, that is £(ko,loq) = 0. Assume in addition 
the generic conditions Im£ 2 7^ 0, ^(kq,loq) 7^ 0, ^(kq,ioq) 7^ 0, ^(kq,loq) 7^ hold. 
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Figure 11: Schematic of the functions ui a and Lot for two cases Mo = 2, Mi = 1, k = 1, 71 = 7. The solid dot 
represents where the guided modes exist. Left: At Kq — for 70 = 7o- Right: At k = =f0.0616 for 70 = 1. 



1. There exist intervals I about kq and V about ujq and analytic real-valued functions uj a ,ujb '■ 
I — > V such that a(n,uj a (K)) = 0, 6(k, w?,(k)) = and u a (Ko) = ujb(^o) = ^o- Thus uj a (n), 

f or K G A{ K o} describe real frequencies for which the transmission T reaches presicely 
100% (peak) and 0% (dip), respectively. 

2. Either uj a {K) > oj^h), for all k G /\{ko}; which means the peak in the transmission comes to 
the right of the dip, or ui a {n) < uJb{^), K G I\{kq}, which implies the peak in the transmission 
comes to the left of the dip. 



Proof: According to (3.47) the zero-sets for a(n, oj) and b(K,uj) are defined by 



det 






-To 



-(7o + 7i) 



V -71 



UJ 



UJ 



2k 
M 
2k 
Mi 



71 

2k cos(7tk) 



v / MoMi" 
2k cos(ttk) 





-2i sin(27r0i; 
-7o 
7i 



-2 7 o 


2uJ -Mo~ 
Ak cos(7tk) 

yJMoh'h 



\ 



0. 



(5.91) 



det 



UJ 



\ 



UJ 



71 

2k 

M 
2k 

M-i 



70 

2k cos(7tk) 

VM0M1 
2k cos(ttk) 



4i sin(27r6»i) 
-7o 
7i 



-2 7 o 

4fc cos(7tk) 
y/M Mi 



0, 



(5.92) 



respectively, which arc real-valued functions of (k.,u>) in region I of the real (k,u>) plane, with 
sin(27r#i) = iyj (2 — | + cos(7tk)) 2 — 1. Both of these conditions are satisfied at (k,q,ojq). Because 

of the conditions ^|(k ,w ) 7^ an d J^I(/t ,wo) 7^ ^' P ar ^ W °f t ne theorem follows from the implicit 
function theorem. These functions have expansions with real coefficients, 



=u -£i(k- kq) - r 2 (n - k Y 

W 6 (k) =0J -£ 1 (k- Kq) -t 2 (K- Kq) 2 



(a = 0) 
(6 = 0) 



(5.93) 
(5.94) 



in which the coefficients of the linear terms are equal by Theorem [15] 

In |16l Theorem 20(4)], it is proved that the assumption Im^2 7^ implies that r 2 7^ t 2 , from 
which part (2) follows. I 

The analysis of the transmission anomaly relies on the following conditions: 



(k,o;,7o)| 2 = \a(K,uj,j )\ 2 + |6(/c, w, 7o)| 2 for k,uj,j G 
if ^(k,w,7o) = for k G R, then Im(cj) < 0. 



(5.95) 
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£(«5,Wo.7o) = 0, a(«o»w5 > 7o) = 0, &(kS,o>S,7o) = (5.96) 

where («g = 0> w cS)7o) ^ ^ 3 * s * ne bifurcation point. 
The following conditions hold generically: 

^(«S,w5, 75)^0, ^K,^,7o*)^0, £( K *, w *, 7 *)/0. (5.97) 

The curves a(re, w, 70) = and b(n, to, 70) = for real values of k near the bifurcation point describe 
frequencies w , cj& of the reflected and transmitted coefficients, respectively, which correspond peaks 
and dips of the transmission. 

Theorem 18. Suppose that for the period N = 2 system with fixed real values of Mq, M\, h\ = 
hi = k, and 71, there exists a unique triple (kq = 0,Wq,7o) G ^ 3 with (wqj7o) ^ n ^ e regime of 
one propagating harmonic, such that £(k,q,Uq,^/q) = 0. Let £{k,uj,^q) = £i(k, cj, 70) + i£2(ft, w, 70); 
where C\ = Re(^), £2 = Im(£) and £1, £2 are real- analytic functions of the real triple (k,w,7o). 
Assume (5.97) ZioZd and 

det f ) ^0. (5.98) 

T/ien there are intervals I about k* q , J about 7 q, and V about ojq and smooth real-valued functions 
u a , coi, : I x J — >■ V , g : I — )> J W : I V such that 

a(K,w a (K,7o),7o) = and &(«, a; 6 («, 70), 7o) = for all (k,7 ) G I x J; 
£(kq, W(ko),5(ko)) = /or aZZ k G /, 
ff(«o) = 7o anrf w ( K *o) = w o- 

Let its make the generic assumption that </'(kq) < (resp. </'(avq) > 0) and that uo a {K,^Q) < 
w &( k >7q) / or some K E I (an analogous conclusion holds for ">"). 
The system undergoes a bifurcation at 70 = 7q •' 

1. For 70 = 7qj there is a unique kq G I such that g(no) = 70, namely, kq = Kq = 0. Moreover, 
W(k* ) =uj a (K* ,%) = w fc (K*,7*) = u*; and w a (K, 7 *) < uj b {n,%) for all k G I\{k*}. 

2. For each 70 G J with 70 < 7g (resp. 70 > 7g^, i/tere exists exactly one kq > in I such 
that 70 = g(±Ko), w a (±K ,7o) = ^£>(±^o,7o) = W(±k ); and w a (K,7 ) < w fe (K,7 ) for all 

K G J\{-K ,K }. 

5. For eac/i 70 G J u>z£/i 70 > 7q (Vesp. 70 < Jq), there exists no kq in I such that 70 = g{no); 
and uj a (n, 70) < ^(k, 70) /or all k G I. 

Proof: The existence of the stated intervals and real-analytic functions w a &> 9> an d is a con- 
sequence of the implicit function theorem. Because of the symmetry of £, a, and b in k, both g 
and are also symmetric. This and the nonvanishing of the second derivative ^"(kq) give rise to 
the three cases depending on 70. The equality u a (no, g(no)) = ^(kq, g(no)) = W(kq) comes from 



the conservation of energy relation \a\ 2 + |6| 2 = \£\ 2 for real (re, uj). Lemma 17 guarantees that 
Wa(^>7o) — w b( K )7o) for all re G I or w a (re,7o) > C0(,(k, 70) for all re G J. I 

The transmission coefficient near the bifurcation point depends delicately on the three analytic 
parameters 00, re, and 70. To obtain an asymptotic formula for the transmission anomaly near 
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the bifurcation, we use, in place of 70, the wavenumber kq of a guided mode at an isolated pair 



(ko,wo) in real (k, w)-space; indeed, Theorem 18 tells us that 70 is an analytic function of kq. We 
do a complex-analytic perturbation analysis in the variables (k = k — k,q,uj = uj — ujo,kq) about 
(0, 0, 0), keeping in mind that ujq = W(k,q) depends analytically on kq. The Weierstrafi preparation 
theorem for analytic functions of three variables provides the following expansions for I, a, and b 
near (/£, uj, kq) = (0, 0, 0): 



£ = e m [u + £ lfiKQ + t Ql n + £ hl K K + ...+ iij^W + 

x [Ao + Xiuj + A 2 ko + + 
w + riflKo + ro,i« + ri f i«o« + . . . + t^kqW + 

X [p + pi5 + p 2 K + /9 3 K + 
t ,lK + tl,l«0« + • • • + tijKQ l W + 

x [to + + T2K0 + r 3 /5 + 



,#2 



e^ 3 [a5 + tin«o 



Taking into account the symmetry of these functions in k and kq, we obtain 



£ = e^i [ w + 4 )lKoK + 4 + ...][! + l x uj + L 2 k k + + . . .] 



a 
b -- 



[1 + Piw + P 2 kok + P3K 2 + 



e^ 3 ro[w + ti,i« K + *o,2« +...][! + Qiw + + Q3K + • • •] 



(5.99) 
(5.100) 
(5.101) 



(5.102) 
(5.103) 
(5.104) 



Inserting these expressions into the law of conservation of energy for real (k,uj,kq) and matching 
like terms yields relations among the coefficients; for example, 



(uj 2 term) 
(ujkok term) 
(ko 2 k 2 term) 
(ujk 2 term) 
(k 4 term) 



1 



Po 2 + r 2 



Re(4,i) 
14,1 1 2 = 
Re(4, 2 ) 



0,2 1 



= p5Re(ri j i) + 7ffRe(ti,i) 
Po 2 |ri,i| 2 + r 2 |ti i i| 2 , 
= plRe(r 0j2 ) + r 2 Re(t ,2) 
Po\r ,2\ 2 + r 2 |t ,2| 2 - 



(5.105) 



Figs. 12 and 13 show the transmission coefficient at and after the bifurcation, both by direct 
calculation as well as using the above expansions in the expression |6|/|^| = H/\/l a l 2 + l^l 2 w hh 
appropriate choices of coefficients up to quadratic order. 



5.3 Resonant Enhancement 

The transmission anomalies that we have analyzed are accompanied by the perhaps more elemen- 
tary phenomenon of resonant enhancement of the amplitude of the field in the waveguide. This 
enhancement will be manifest in the resonant component of the field X , which is the first term of 



equation (5.79), 



X = caCe 1 + c£C(0,M- l F 2 ) 



(5.106) 



Thus a good measurement of amplitude enhancement is the ratio \a/£\. Let a in the vicinity of 
(kq,ujo) have the expansion 

a = /3 + /3 1 K + /3 2 w + --- (5.107) 



Theorem 19. In the expansion (5.107), /3q = 
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Figure 12: At bifurcation: The transmission coefficient for M = 2, Mi = 1, ki = k 2 = 1, 70 = 7o = 1.029633513, 
71 = 7 near the parameters of a guided mode, kq = 0, ui^ — 0.9778859328. Upper: Direct calculation of the 
transmission coefficient. Lower: The asymptotic formula for 



Proof: By Theorem 10, at the pair (kq, ojq) there is a solution to the scattering problem BX = F, 
with F = C(a, F2). Using £(kq,ujq) = 0, this equation becomes 



(5.108) 



It follows that a(Ko,LOo) = 0. I 



Using the forms for a and (5.81) for I , we obtain 

PlK + (3 2 u H 



n 



(to + £ 1 k + £ 2 k 2 H ) V e * Pl 



(5.109) 



When k is small, the magnitude of the denominator in a/£ is minimized to order 0(k 2 ) when 
see the response to an incident plane wave at this optimal frequency, put 



ui = — l\K — Re(^2)K 2 , or u> = ojq — £{k — Re^)^ 2 , 
and obtain for the amplitude enhancement A 

Pi - hfa ~ f3 2 Re(£ 2 )7i + • ■ 



A 



a 


1 


1 


K 



ilm(£ 2 ) + 

so that A has an asymptotic expansion of the form 

di 



A 



~ — + d 2 + 

K 



(2 



iK-Yle{£ 2 )K A , K-tO) 



(5.110) 



(5.111) 



(5.112) 
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Figure 13: After bifurcation: The transmission coefficient for Mo = 2, M\ = 1, ki = fea = 1, 7o — 7o — 
0.0001 = 1.029533513, 71 = 7 near k — and near the parameters of the guided mode, ko = 0.003564296929, 
cjo = 0.9778903229. Upper: Direct calculation of the transmission coefficient. Lower: The asymptotic formula for 
\b\/\i\. 
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Figure 14: The dots represent numerically calculated magnitudes of the field in the waveguide produced by an 
incident plane wave of amplitude 1 at various values of k and at uj = uj — l-Ji — Re(^2)£5 2 for two cases of period 
N = 2. Left: k* = 0, with M = 2, Mi = 1, k = fci = 1, 70 = 1.029633513, 71 = 7. Right: k = 0.0616, with 
M = 2, Mi = 1, k = k 1 = 1, 70 = 1, 71 = 7. 
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Fig. 14 shows numerical calculations that confirm the behavior of the amplitude of the waveg- 
uide at to = —iik — Re(^2)'5 2 for N = 2. The magnitude of the field is calculated using the 
expression y^col 2 + |ci| 2 . 

6 Appendix: Difference Operators 

We use the following notation: 

V = {Vn}, W = {w n }, VW = {v n W n } , 
{v + ) n = V n+ i, {v~) n = V n -X, 

{d x v) n = {v x )n = V n+ i - V n , {d x v) n = {v x ) n = V n - U n _l- 

One can compute the discrete product rule and the fundamental theorem as well as a summation- 
by-parts formula: 

(vw) x = v x w + + vw x = v x w + v + w x , (6.113) 
(vw)x = v x w~ + vw x = v x w + v~w x , (6.114) 

m 2 m 2 — 1 

^7712 %i — 



^ (Vx) m = ^2 (^) m ' (6.115) 
m=mi+l m=m\ 
m 2 m 2 

^ {d x v) m w m + ^ v^d x w m = v m2 w m2 - v mi w mi . (6.116) 

m=m\ + l m—mi+1 

For the two-dimensional discrete calculus, we use the notation 

F = {F mn } = {(F mn , F mn )} , 

(F )mn = (-^m— X,m F m,n— l)i 

v- = (d x ,dy), v + = (d x ,d y ). 

One has the discrete product rule 

V_ • (vF) = v\7- ■ F + V_u • F . (6.117) 
The discrete divergence theorem for a rectangular region [mi, m.2] x [^1,^2] is 

n 2 m 2 n 2 m 2 

£ £ ( V -- F U= £ ( i? m 2 n-^ 1 J+ E ( F L 2 -Fl ni ). (6.118) 

n=ni+l m=mi + l n=ni+l m=mi+l 



Now put F = vV+u in (6.118) and expand V_ • (vV + u) using (6.117) to obtain the two-dimensional 



summation-by-parts identity 



TI2 m 2 n 2 

n=ni+im=mi+l n=n\+l 



ni2 

+ E (Wmnj-KWJ (6.119) 
m=mi+l 

n 2 m 2 

n=n\+l m=mi+l 
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Here, one uses the identities V_ • V + u = An and (V+it) = V_u. 
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